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Abstract 

The ®B solar neutrino flux that has recently been measured by Super- 
Kamiokande is consistent with the ®^Ar production rate in ^^Cl at Homes- 
take. The gallium solar neutrino experiments, GALLEX and SAGE, con¬ 
tinue to observe ^^Ge production rates in ^^Ga that are consistent with the 
minimal signal expected from the solar luminosity. The observed ®B solar 
neutrino flux is in good agreement with that predicted by the standard so¬ 
lar model of Dar and Shaviv with nuclear reaction rates that are supported 
by recent measurements of nuclear fusion cross sections at low energies. 
But, the signals measured by Super-Kamiokande, SAGE and GALLEX 
leave no room for the contribution from the expected ^Be solar neutrino 
flux. This apparent suppression of the ^Be solar neutrino flux can be 
explained by neutrino oscillations and the Mikheyev-Smirnov-Wolfenstein 
effect, although neither a flavor change, nor a terrestrial variation, nor a 
spectral distortion of the ®B solar neutrino flux has been observed. De¬ 
tailed helioseismology data from SOHO and GONG confirm the standard 
solar model description of the solar core, but helioseismology is insensi¬ 
tive to the fate of ^Be in the sun and the production rates of ^Be, CNO 
and ®B neutrinos. Thus, helioseismology and the solar neutrino problem 
do not provide conclusive evidence for neutrino properties beyond the 
standard electroweak model. The solar neutrino problem may still be an 
astrophysical problem. The deviations of the experimental results from 
those predicted by the standard solar models may reflect the approximate 
nature of our knowledge of nuclear reaction rates and radiation trans¬ 
port in dense stellar plasmas and the approximate nature of solar models. 
Only future observations of spectral distortions, or terrestrial modulation 
or flavor change of solar neutrinos in solar neutrino experiments, such as 
Super-Kamiokande, SNO, Borexino and HELLAZ will be able to establish 
that neutrino properties beyond the minimal standard electroweak model 
are responsible for the solar neutrino problem. 
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1 Introduction 


The sun is a typical main sequence star that generates its energy via fusion of 
hydrogen into helium through the pp and CNO nuclear reaction chains (Fig. 1). 
Due to the conservation laws of baryon number, electric charge, lepton flavor and 
energy, the total solar neutrino flux is fixed, practically, by the solar luminosity 
(e.g., Dar and Nussinov 1991). The neutrino spectrum is essentially a sum 
of standard beta decay spectra from the /3-decays 2p ^ De+t'e, ® 2ae+^e, 

and and “lines” from the electron captures e^Be^ 

j^Jbi and pep ^ Dr's- To a good approximation, they are independent of the 
conditions in the sun. However, the relative contributions of the various solar 
neutrino sources depend on the chemical composition, temperature and density 
distributions near the center of the sun. These are usually estimated from 
standard solar models (SSM). They can be tested also by helioseismology, which, 
however, is not sensitive to the exact abundance or fate of trace elements (e.g., 
^Be) in the sun. 

Solar neutrinos, have been detected on Earth in roughly their expected num¬ 
bers, in five underground solar neutrino experiments, the Chlorine solar neu¬ 
trino experiment at Homestake, South Dakota, USA, the Water Cherenkov Ex¬ 
periment, Kamiokande, at Kamioka Japan, the Soviet-American Gallium Ex¬ 
periment, SAGE, at the Baksan, Russia, the European Gallium Experiment, 
GALLEX, at Gran Sasso, Italy and the large water Cherenkov Experiment, 
Super-Kamiokande, at Kamioka Japan. These experiments have confirmed that 
the sun is powered by fusion of hydrogen into helium. This milestone achieve¬ 
ment in physics, however, has been overshadowed by the fact that the combined 
results from the solar neutrino experiments seem to suggest that the solar neu¬ 
trino flux differs significantly from that expected from the standard solar mod¬ 
els. This discrepancy has become known as the solar neutrino problem. Many 
authors have argued that the solar neutrino problem provides conclusive evi¬ 
dence for neutrino properties beyond the minimal standard electroweak model. 
However, conclusive evidence for new electroweak physics from solar neutrino 
observations can be provided only by detecting at least one of the following 
signals: 

1. Solar neutrinos with flavors other than Vf.. 

2. Spectral distortion of the fundamental /3-decay spectra. 

3. Terrestrial modulations of the solar neutrino flux. 

4. A clear violation of the luminosity sum rule. 

So far, no such conclusive evidence has been provided by the solar neutrino ex¬ 
periments. Therefore, the solar neutrino problem does not provide solid evidence 
for neutrino properties beyond the standard electroweak model and standard 
physics solutions to the solar neutrino problem are not ruled out. Moreover, 
a closer look at the sun through helioseismology. X-ray and UV observations 
shows that the sun is a bewildering turmoil of complex phenomena. It shows 
unexpected features and behavior at any scale. It has a strange complex in- 
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ternal rotation, unexplained magnetic activity with unexplained 11 year cycle, 
unexpected anomalies in its surface elemental abundances, unexplained explo¬ 
sions in its atmosphere and unexplained mechanism that heats its million degree 
corona and accelerates the solar wind. Perhaps the surface of the sun is com¬ 
plex because we can see it and the center of the sun is not only because we 
cannot? Perhaps the SSM which has been improved continuously over the past 
three decades but which still uses simple plasma physics and assumes an exact 
spherical symmetry, no mass loss or mass accretion, no angular momentum loss 
or gain, no differential rotation, zero magnetic field through the entire solar 
evolution, is a too simplistic picture and does not provide a sufficiently accurate 
description of the core of the sun. 

In this paper we summarize the experimental results from the various solar 
neutrino experiments, we discuss shortly the solar neutrino problem in a “model 
independent way” and we compare the experimental results with updated stan¬ 
dard solar model calculations. We conclude that (a) there is no conclusive 
evidence for a ®B solar neutrino problem, (b) the ®B solar neutrino flux as 
measured by Super-Kamiokande is in good agreement with that predicted by 
the standard solar model (of Dar and Shaviv) with nuclear reaction rates that 
are supported by recent measurements of nuclear reaction rates at low ener¬ 
gies, (c) the suppression of "^Be solar neutrinos which is suggested by both the 
chlorine and gallium experiments can be due to neutrino oscillations and the 
Mikheyev-Smirnov-Wolfenstein effect, although neither a day-night effect, nor 
a spectral distortion of the ®B solar neutrino flux, nor terrestrial modulation 
of the flux, has been observed, (d) a deficit of ^Be solar neutrinos, if there is 
one, may still be explained by standard physics and/or astrophysics (e) only 
future observations of spectral distortions or flavor change of solar neutrinos 
in solar neutrino experiments, such as Super-Kamiokande, SNO, Borexino and 
HELLAZ may establish that neutrino properties beyond the minimal standard 
electroweak model are responsible for the solar neutrino problem. 

2 Solar Neutrino Observations 

2.1 The Chlorine Detector 

The radiochemical chlorine detector (Davis 1966, Cleveland et al. 1998) which 
contains 615 tons of tetrachloroethylene, C2CI4, measures the production rate of 
^^Ar by solar neutrinos through the reaction I'g -|-^^ Cl Ar + e~ which has 
a threshold energy of 814 keV. The detector is located at the Homestake gold 
mine at a depth of 1480 meters underground (~ 4200 meters water equivalent) 
in Lead, South Dakota, USA. After an exposure time of 1 to 3 months the 
^^Ar atoms are extracted from the target liquid and counted by observing their 
electron capture decay back to ^"^Cl (half life T 4/2 = 35 d) in a proportional 
counter. The ^^Ar production rate measured in 108 individual runs (runs 18 to 
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133) from 1970 to 1995 is plotted in Fig. 2. The mean capture rate of solar 
neutrinos was 


< cr(/)^e >ci= 2.56 ± 0.16(stat) ± 0.14(syst) SNU (1) 

where 1 SNU=1 neutrino capture per second in 10^® target atoms. 

2.2 The Light Water Cherenkov Detectors 

Kamiokande II, the first imaging light water solar neutrino detector, (later up¬ 
graded to Kamiokande III) was located in the Kamioka mine in Japan (Fukuda 
et al. 1996). It measured the Cherenkov light emitted by electron recoils pro¬ 
duced by elastic scattering of solar neutrinos from electrons in the inner 680 tons 
of a large tank filled with a total of 2180 tons of light water. For radioactive 
background reduction the threshold had to be set to a rather high electron recoil 
energy (7 MeV). Therefore, the detector was sensitive only to the upper end of 
the ®B solar neutrino spectrum. During cumulative live time of 2079 days the 
detector yielded an average flux of 

<('i/e(*B) = 2.80 ± 0.19(stat) ± 0.33(syst) x 10® cm“^s“^, (2) 

as shown in Fig. 3. Within its experimental sensitivities KII-I-KIII has not 
detected temporal variation or spectral distortion of the the ®B spectrum. 
Super-Kamiokande (SK) is a 500,00 tons imaging water Cherenkov detector 
whose inner 22,500 tons are used for the solar neutrino measurements (Fukuda 
et al 1998). The current energy threshold is 6.5 MeV. It started its operation on 
April 1 1996. Fig. 4 shows the distribution of events as function of the cosine of 
the angle of electrons recoiling from neutrino-electron scattering relative to the 
direction from the sun during live time of 504 days. The solid line is the best 
fitted histogram due to the ®B solar neutrino flux of (Suzuki et al. 1998) 

<?ii/e(®B) = 2.44to;Q5(stat)tQ;g7(syst) x 10® cm“^s“^ (3) 

and a constant background. The observed energy spectrum of the recoiling 
electrons is consistent with that expected from elastic scattering of ®B solar 
neutrinos. In spite of the large number of events the statistics are not large 
enough yet to show the expected « 6.5% periodical variation in the flux due to 
the periodical variation in the distance of Earth from the sun, as can be seen 
from Fig. 5. The SK data also does not show any dependence on the path 
length of solar neutrino in Earth, as can be seen from Fig. 6. The SK limit on 
a Day-Night asymmetry is 

D - N 

A = = -0.023 ± 0.020(stat) ± 0.014(syst). (4) 
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2.3 The Gallium Detectors 


The radiochemical gallium detectors measure the production rate of ^^Ge by so¬ 
lar neutrinos through the reaction t'e Ga Ge -I- e“. The energy thresh¬ 
old, 232.2 keV, is well below the maximum energy of the pp neutrinos, 420 keV. 
After exposure time of a couple of weeks the ^^Ge atoms are extracted from 
the target liquid and counted by detecting the X-ray emission from their elec¬ 
tron capture decay back to "^^Ga (half life T 1/2 = 11-43 days) in a proportional 
counter. 

The Gallex detector (Hampel et al. 1996; Kirsten 1998) is located in the 
Gran Sasso underground laboratory in Italy. It contains 30.3 tons of gallium 
in GaGls — HCl solution. The neutrino produced ^^Ge atoms form the volatile 
compound GeCh which at the end of the exposure is swept out of the solution 
by a gas stream and converted into GeH 4 . The produced ^^Ge are counted by 
detecting their radioactive decay inside a proportional counter. The combined 
results of 65 individual GALLEX runs corresponding to data taking periods 
GALLEX I,II,III and IV as shown in Eig. 8 gave a production rate of 

< >Ga= 76.4± 6.3(stat)t49(syst) SNU. (5) 

GALLEX has also conducted two ^^Gr neutrino source experiments to test the 
overwhole performance of the detector. The ratio between the measured ^^Ge 
production rate due to the ®^Cr source and the expected rate from the known 
source strength was 1.00 ±0.10 and 0.83 ±0.10, in the two experiments, respec¬ 
tively. 

The Soviet American Gallium Experiment (SAGE) is located in the Bak- 
san Neutrino Observatory in the northern Caucasus of Russia at a shielding 
depth of 4715 meter water equivalent. The detector uses 55 tons of metallic 
gallium. After the metal is converted to solution the produced ^^Ge atoms are 
removed from the Gallium and counted in a low background proportional coun¬ 
ters by a procedure similar to that used by GALLEX. The combined results of 65 
individual runs of SAGE corresponding to data taking periods SAGE I,II,III,IV 
as shown in Fig. 8 gave a production rate of (Abdurashitov et al 1996, Gavrin 
1998) 

< <y4>v<, >Ga= 70 ± 6.3(stat)t4;g(syst) SNU. (6) 

SAGE has also conducted a ®^Cr neutrino source experiment. The ratio between 
the measured "^^Ge production rate due to the ^^Cr source and the rate expected 
from the source strength was 0.95 ±0.11. 

3 The Luminosity Sum Rule 

Due to conservation of baryon number, electric charge, lepton flavor and energy, 
the net reaction in the sun can be written as 

2e- ± 4p ^4 jje ± 26.732 MeV. (7) 
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Thus, the generation oi Q = 26.732 MeV in the sun is accompanied by the 
production of two Ve^s. If the sun is approximately in a steady state with a 
nuclear energy production rate that equals its luminosity, then the total solar 
neutrino flux at Earth is (Dar and Nussinov 1991), 

2 S' 

where Sq = Lq/AttDq « 1.367 kW is the measured “solar constant” 
which yields a solar luminosity Lq « « 3.846 x 10^^ erg s~^ for an 

average distance Dq « 1.496 x 10^^ cm of Earth from the sun, and = 
/(j)^^ is the mean energy of solar neutrinos which has been ap¬ 
proximated by (pp) « 0.214 MeV, the mean energy of the pp solar neutrinos 
that dominate the solar neutrino flux. Eq. (2) can also be rewritten as a sum 
rule, 

S,(Q/2-E,,(i))</),^(,)«5e. (9) 

The summation extends over all the neutrino producing reactions with = 
0.265, 1.442, 0.814, 6.710, 0.707, and 0.997 MeV for the pp, pep, "^Be, ®B, 
and neutrinos, respectively. If the small pep flux, which is proportional to 
the pp flux, is included in the pp flux, and if the very small ®B neutrino flux is 
neglected then the solar luminosity sum rule can be rewritten as 

0.9800(/)^3(pp) -f 0.939(?ii,,(Be) -b 0.936(?ii,,(NO) « 6.377, (10) 

where are in units of 10^° cm“^s“^. 


4 Helioseismological Constraints 

Accurate ground based (e.g., Hill et al. 1996 and references therein) and space 
based measurements aboard SOHO (e.g., Turck-Chieze et al. 1997 and refer¬ 
ences therein) of solar photospheric oscillation frequencies provided detailed in¬ 
formation about the structure of the solar interior (e.g., Christensen-Dalsgaard 
1996). In particular the base of the convection zone has been determined to 
be at Rcz « (0.713 ± O.OO3)R0 (Basu and Antia 1997) and the photospheric he¬ 
lium abundance has been inferred to be Yg = 0.249 ± 0.003 (Basu and Antia 
1995). Helioseismology is generally in good agreement with the standard solar 
models (see e.g. Christensen-Dalsgaard 1996). However, there are systematic 
deviations between the helioseismology determination of the sound speed in the 
sun and that predicted by the SSM (see, e.g., Eig. 10.) which are similar to all 
SSMs and whose origin is not clear. Moreover, helioseismology is sensitive only 
to the average local properties in the sun (temperature, density, average molec¬ 
ular weight) but not to the rates of the rare nuclear reactions in the sun which 
produce the pep, ^Be, ®B and NO solar neutrinos. Therefore, helioseismology 
should not be used to argue that the SSM predict correctly these solar neutrino 
fluxes. 
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5 Model Independent considerations 

Counting rates in solar neutrino experiments are formally given by 

= / idn,JdE)a,a{E)dE (11) 

JEo 

where Na is the number of “active” atoms in the detector, a^a{E) is their 
cross section for neutrinos with energy E, Eq is the threshold energy of the 
detector, dn,^JdE is the energy spectrum and 4>ue{i) is the total flux of neutri¬ 
nos from reaction i in the sun. Both, dn,yJdE and cr^a follow from the stan¬ 
dard electroweak theory and are essentially independent of the sun: dn^^/dE 
is practically the standard /3-decay spectrum for the /3-decays 2p^ De^Ve, 
®B—> 2ae“'':^e5 and and is a J-function for the 

electron captures e“^Be^ ^Li and pep Due.) Thus conclusive evidence 
for new electroweak physics can be provided only by detecting at least one of 
the following signals: (1) Spectral distortions of the /3-decay spectra of solar 
neutrinos. (2) Solar neutrino flavors other than v^- (3) Terrestrial modulations 
of the solar neutrino flux. (4) A clear violation of the luminosity sum rule. 

So far, no such conclusive evidence has been provided by the solar neutrino 
experiments. 

1. Spectral Distortion. At present only Super-Kamiokande can measure 
the spectrum of solar neutrinos (above 6.5 MeV). Within their statistics and 
systematic uncertainties the energy distribution of the detected electrons which 
are scattered by solar neutrinos is consistent with that expected from an undis¬ 
torted spectrum of ®B neutrinos. This can be seen from Fig. 7. A “hint” for a 
spectral distortion may exist in the SK data, but it depends strongly on events 
beyond the kinematical limit which are attributed to the detector energy res¬ 
olution. Also the neutrino spectrum near the “end point” is not known very 
well because the ®B decays into a virtual short lived 2 q; state which has a large 
energy spread. 

Super-Kamiokande, which has been running since April 1, 1996, will finally 
have a much larger statistics and perhaps a lower threshold energy. These, 
perhaps, will be able to provide more conclusive evidence. 

2. Neutrino Flavor Change. The radiochemical experiments are blind to 
neutrino flavors other than that of VeS. SK is sensitive also to t'^’s and VtS 
but it cannot distinguish between the neutrino flavors at solar energies. Only 
future experiments, such as SNO, will be able to obtain information on the 
flavor content of the solar neutrino flux. 

3. Terrestrial Modulations. In spite of the large number of events collected 
by SK the statistics are not large enough yet to show the expected « 6.5% 
periodical variation in the flux due to eccentricity of the Earth’s orbit around 
the sun, as can be seen from Fig. 5. The SK data also show no dependence 
on the path length of solar neutrino in Earth, as can be seen from Fig. 6, no 
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day-night effect and no winter-summer difference. In particular, its limit on the 
Day-Night asymmetry is A = (D — N)/(D -|- N) = —0.023 ± 0.025. 

4. Violation of the Luminosity Sum Rule. A clear violation of the so¬ 
lar luminosity sum rule could prove that lepton flavor is not conserved. The 
“minimal” expected signal in gallium which follows from the luminosity sum 
rule is obtained by assuming that all the solar neutrinos are pp neutrinos. 

If the mean cross section for the capture of the pp neutrinos in gallium is 
a ~ (1.17 ± 0.03) X 10"^® it yields a minimal signal of 76 ± 2 SNU in 

^^Ga. In fact, the solar neutrino flux, which is observed by SK contributes 
additional < a(j){^B) >Ga= 5.85 ± 1.75 SNU to the ^^Ge production rate and 
the minimal expected signal in ^^Ga is 82 ± 3 SNU. It is consistent, within the 
experimental and theoretical uncertainties, with the 76.4 ± 8 SNU production 
rate of ^^Ge by solar neutrinos in gallium which was measured by GALLEX. 
The gallium experiments, however, appear to leave no extra room for the con¬ 
tribution from ^Be solar neutrinos. 

A further indication that the ^Be and perhaps also the GNO neutrinos are 
missing is provided by the Cl experiment. The expected production rate of ^^Ar 
in ^^Cl by the ®B solar neutrino flux that was measured by Super-Kamiokande is 
< >ci= 2.68 ± 0.15 SNU. It is consistent with < >ci= 2.56 ± 0.21 SNU, 
the total ^^Ar production rate in ^"^Cl, as measured at Homestake but leaves no 
room for the contribution from the ^Be solar neutrino flux. 

The SK measurement (3) and the solar luminosity sum rule (10) two obser¬ 
vational constraints on the solar neutrino fluxes. Two additional constraints are 
provided by the gallium and chlorine experiments. Assuming that the neutrino 
capture cross sections are well represented by their theoretical estimates (e.g., 
Bahcall 1989; 1998 and references therein) one can write them approximately 
as 

Ga : (11.7 ± 0.3)^;.^^ (pp)-f (71.7 ±5.0)())i,,(^Be) -k (2.40 ±0.78) X 10Vi^,(®B) 

± (87± 12)(()^,(NO) = 77.5U7.8 (12) 

Cl : (1.11 ± 0.04) X lOV.e (®B) ± (2.4±)(/),, (^Be) 

± (16±)0,,(pep) ± (4.2±)(/),,(NO) = 2.56 ± 0.25 (13) 

If the theoretical estimates of the cross sections for nuclear capture of solar 
neutrinos represent well their true values then the only physical solution of the 
four constraints is 

(pp + P^P) ~ 6-5 X 10^°cm“^s“^, (®B) « 2.44±0.11 x 10®cm“^s“^, 

("Be) « C“("Be), (NO) «^^(NO). (14) 

The confidence level of this solution cannot be quantified in a reliable way be¬ 
cause of the unknown origin and nature of all the systematic errors in both the 



theoretical cross sections and the experimental results. The only reliable conclu¬ 
sion is that (/)i/^(^Be) and (j),^^(NO) appear to be strongly suppressed compared 
with their standard solar model estimates. 


6 Is There a Solar Neutrino Problem? 

Table I presents a comparison between the solar neutrino observations and the 
SSM predictions of Bahcall and Pinsonneault 1995 (BP95), Bahcall, Basu and 
Pinsonneault 1998 (BSP98), Brun, Turck-Chieze and Morel 1998 (BTM98) and 
Dar and Shaviv 1996 (DS96). Although BP95 and BSP98 predict a ®B solar 
neutrino flux that is approximately 2.2 and 2.7 times, respectively, larger than 
that observed by SK, DS96 predict a flux consistent with that observed by SK. 
The differences between BP95, BSP98 and DS96 are summarized in Table II 
(for details see Dar and Shaviv 1996). The difference between the predicted 
®B flux are mainly due to the use of different nuclear reaction rates in DS96, 
differences in the calculated effects of diffusion and differences in the initial 
solar chemical composition assumed in the two calculations. They reduce the 
predicted ®B flux relative to those in BP95 (BSP98) by approximate factors of 
0.55 (0.70), 0.81, and 0.95, respectively. The remaining differences are mainly 
due to inclusion of partial ionization effects, premain sequence evolution and 
deviations from complete nuclear equilibrium in DS96 which were neglected in 
BP95 and BSP98, and due to different numerical methods, fine zoning and time 
steps used in the two calculations: 

6.1 Initial Chemical Composition 

The initial chemical composition of the sun influences significantly the solar 
evolution and consequently the present density, chemical composition and tem¬ 
perature in the solar core that determine the solar neutrino fluxes. In particular, 
the radiative opacities, which determine the temperature gradient in the solar 
interior, are very sensitive to the abundance of heavy elements which are not 
completely ionized in the sun. Since the initial chemical composition of the 
sun is unknown, one must infer it indirectly, e.g., from the chemical composi¬ 
tion of the solar photosphere, the primitive early solar meteorites and the local 
interstellar medium (ISM). 

The solar photospheric abundances have changed only slightly during the 
solar evolution by gravitational settling, diffusion and turbulent mixing in the 
convective layer and by cosmic ray interactions with the surface material. There¬ 
fore, in principle, one can adjust the initial chemical composition of the sun in 
the SSM to yield its measured photospheric composition. Unfortunately, the 
photospheric abundances of most elements are still not known to sufficient ac¬ 
curacy and there is no direct spectroscopic information on the photospheric 
abundance of "'He. Consequently, the initial mass fraction Yi of ^He in the sun 
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has been treated in the SSM as an adjustable parameter. Recently, however, 
the photospheric mass fraction of ^He has been inferred from helioseismologi- 
cal measurements of the sound speed in the convective solar layer. The best 
estimated value is now Yg = 0.249 ± 0.003 (Basu and Antia 1997). Since the 
photospheric mass ratio of metals ((A > 4) to hydrogen is (Grevesse, Noels and 
Sauval 1996) Z/X = 0.0244, one obtains that Xg = 0.733 and Zg = 0.0178. It is, 
however, important to note that the metallicity in the present local ISM, essen¬ 
tially measured from analyses of the Orion nebula and of nearby B stars (Gies 
and Lambert 1992; Cunha and Lambert 1992; Wilson and Rood 1994; Mathis 
1996) is lower than the value obtained from the photospheric abundances. This 
is in contradiction with galactic chemical evolution models which predict an in¬ 
crease with time of metallicity in the ISM. Moreover, gravitational settling and 
diffusion decrease the solar surface metallicity with time. 

The meteoritic elemental abundances are known with much better accuracy. 
Beside the noble gases, and H, G, N and O which were able to form highly 
volatile molecules or compounds and escape condensation, all the other elements 
are believed to have condensated completely in primitive early solar system 
meteorites. Therefore, their initial relative abundances in the sun are expected 
to be well represented by their values in type I carbonaceous chondrites. If 
diffusion and gravitational settling have not changed their ratios significantly, 
then the relative abundances of these elements in meteorites must be similar to 
their photospheric values. Over the past decades there have been many initial 
disagreements between the meteoritic and photospheric abundances. In nearly 
all cases, when the atomic data were steadily improved and the more precise 
measurements were made, the photospheric values approached the meteoritic 
values. The photospheric abundances are now as a rule in very good agreement 
with the meteoritic values if the conversion factor from the solar abundance 
scale Nh = 10^^ to the meteoritic scale. Ns, = 10® is R = log(sol)/met) = 1.56, 
as can be seen from Table III borrowed from Grevesse, Noels and Sauval 1996. 

In our SSM, Yi was left as an adjustable parameter. The initial solar heavy 
metal abundances were assumed to be equal to the meteoritic (Cl chondrites) 
values of Grevesse, Noels and Sauval (1996). The overwhole initial metallicity 
ratio Zj/Xj and the initial H, C, N, O and Ne abundances were adjusted to 
yield their present photospheric values (Zg/Xg = 0.0244 and the values quoted 
in Table III, respectively). Our SSM yields Yg = 0.238 ± 0.05 in agreement with 
the helioseismological estimates. 

The photospheric abundances of ^Li, however is smaller by a factor of nearly 
140 than its meteoritic abundance. The origin of such large difference is still 
not clear. It cannot be explained by nuclear burning during the Hayashi phase 
although significant lithium burning does takes place during this phase. It may 
be explained by rotational mixing (e.g., Richard et al 1996). Although the 
initial solar (meteoritic) abundances of lithium is very small and do not play 
any significant role in solar evolution, its depletion perhaps can provide a clue 
to the real history of the convection zone and the sun. 
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6.2 Nuclear Reaction Rates 

The nuclear reaction rates for most stellar reactions are inferred by extrapo¬ 
lating measurements at higher energies to stellar reaction energies. The cross 
sections at center of mass energies well below the Coulomb barrier are usually 
parametrized as 

a{E) = (15) 

h/ 

where rj = ZiZ 2 e^ jUv is the Sommerfeld parameter, v is the relative velocity 
of the colliding nuclei in the initial state, Zi and Z 2 are their charge numbers 
and E is their center of mass energy. The exponent is an approximate (WKB) 
form for the penetration probability of the Coulomb barrier in the initial state. 
The “astrophysical S factor” is expected to vary only slowly with energy. It 
is usually extracted either from a polynomial fit to experimental data at low 
energies or from theoretical calculations normalized to the experimental data. 

The uncertainties in the nuclear reaction rates at solar conditions are still 
large due to (1) uncertainties in the measured cross sections at laboratory ener¬ 
gies, (2) uncertainties in their extrapolations to lower solar energies, (3) uncer¬ 
tainties in dense plasma effects (screening, correlations, fluctuations and devi¬ 
ations from pure equilibrium distributions) on reaction rates. Rather than av¬ 
eraging measured cross sections that differ by many standard deviations, DS96 
used for their extrapolations only the most recent and consistent measurements 
of the relevant nuclear cross sections. For sub-Coulomb reactions that take 
place when the colliding nuclei are far apart, the Optical Model and the Dis¬ 
torted Wave Born Approximation give a reliable description of their energy 
dependence and can be used for extrapolating the measured sub-Coulomb cross 
sections to solar energies. 

The “astrophysical S factors” which were used in BP95, BSP98 and DS96 
are compared in Table II. The origin of the differences are as follows: 

Sii(O). The value advocated by Adelberger et al. (1998) and used in BSP98 
is based essentially on the updated calculation of Sii(O) by Kamionkowski and 
Bahcall (1994). The authors considered many small corrections. However, they 
ignored the screening effect of the solar plasma electrons in the outgoing chan¬ 
nel. Screening corrections in the nuclear reaction rates were included in BP95 
and BSP98, and in all the published SSM calculations, only in the incoming 
channel. That is justihed for radiative captures like ^He -|-^ He Be -I- 7 and 
p -1-^ Be B -I- 7 because the photon is chargeless. The emitted positron in 
the “beta decay” pp —> D -I- e+ -I- r'e of the two fusing protons sees essentially a 
Deuterium nucleus screened by the plasma electrons in the Debye sphere left by 
the two protons (the mean velocity of the emitted positron is much larger than 
the mean velocity of the electrons in the Debye sphere). If the wave function of 
the ejected positron in the Coulomb held of the Deuteron is calculated from the 
Dirac equation with the Debye screening potential around two fusing protons, 
one obtains that their fusion rate in the sun is enhanced by approximately 1.75% 
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(Dar and N. Shaviv, unpublished). Consequently, the “bare” value of Sii(O) in 
SSM calculations that do not include electron screening in the outgoing channel 
must be increased by 1.75% as in DS96. 

817 ( 0 ). Recent low energy measurement of the cross section for the reaction 
p Be —B + 7 by Hammache et al. (1998) are consistent with the measure¬ 
ments of Vaughn et al. (1970) and of Filippone et al. (1983a,b) which reached 
energies as low as 134 keV. These measurements disagree with the older mea¬ 
surements of Kavanagh (1960), Parker (1966, 1968) and Kavanagh et al. (1969). 

Because of the low binding energy, the radiative capture of p by ^Be takes 
place well outside the range of the nuclear forces. Therefore, at low energies the 
dependence of the cross section on energy is well described by the optical model. 
When applied to the experimental data of Vaughn (1970); Filiponne (1983a,b) 
and Hammache et al. (1998) it yields Sii(O) = 17.8 ± 1.0 eV • b (see also Barker 
1995; Nunes et al 1997). This value is consistent with the indirect measurements, 
Sii(O) = 16.7 ± 3.2 eV • b through Coulomb dissociation of ®B (Motobayashi et 
al 1994) and Sii(O) = 17.6 ± 1 eV • b through proton transfer reactions (Xu et 
al. 1994). The mean value 817 ( 0 ) = 17.5 ± 1.0 eV • b is consistent with the value 
817 ( 0 ) = 17 eV • b used in D896, but it is smaller than both 817 ( 0 ) = 22.4 eV • b 
used in BP95 and 817 ( 0 ) = 19 eV • b that is advocated by Adelberger et al. 1998 
and used in B8P98. 

833 ( 0 ). The cross section for the reaction ^He 4-^ He He -I- 2p has recently 
been measured at the Laboratory for Underground Nuclear Astrophysics (LUNA) 
at energies covering the Gamow peak around Eg = 21.9 keV (Junker et al. 
1998). At such small lab energies screening by atomic electrons enhances con¬ 
siderably the bare nuclear cross section. In the Born-Oppenheimer adiabatic 
approximation, Ue, the gain in kinetic energy by the colliding nuclei, is bounded 
by the change in the total binding energy of the atomic electrons when they oc¬ 
cupy the atomic ground state of the ^He -|-^ He “nuclear molecule”. For ^He^^ 
ions incident on ^He gas target, Ue ~ 240 eV. From eq. (15) one obtains that the 
screening enhancement of the cross section for Ue ^ E is given approximately 

by 


exp^ 


^TrrjUe/E 


(16) 


We stress that one must use a consistent treatment of electron screening 
enhancement of nuclear cross sections in the lab and in the solar plasma. In 
particular, Ue cannot be taken as an adjustable parameter when fitting 8 (E) 
to the lab measurements but must be fixed to its adiabatic value Ug = 240 eV 
in order to be consistent with the “weak screening” prescription that is used in 
the 88 M. For Ue = 240 eV, a best fitted 8 (E) between 20.7 keV and 1080 keV 
to the cross section measurements that are consistent at overlapping energies 
(Dwarkanath and Winkler 1971, Krauss et al. 1987, Greife et al. 1994 and 
Junker et al. 1998) yields the values 833 ( 0 ) = 5.60 MeV • b, 833 ( 0 ) = —4.1b and 
833 ( 0 ) = 4.60 MeV“^ • b. The fit is shown in Fig. 9. Note that if one extracts 
8 (E) directly from the LUNA data at the Gamow peak (Junker et al. 1998) one 
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obtains S 33 (Eg) = 5.75 MeV • b. The value 833 ( 0 ) = 4.99 MeV • b was used in 
BP95, the value S 33 (Eg) = 5.30 MeV • b that was recommended by Adelberger 
et al. 1998 was used in BSP98, while the value 833 ( 0 ) = 5.6 MeV • b was used 
in DS96. 

834 ( 0 ). There are six published measurements of the cross section for the reac¬ 
tion ^He +* He —Be -I- 7 that are based on the detection of the prompt 7 -rays 
(Parker and Kavanagh 1963; Nagatani et al. 1969; Krawinkel et al. 1982; Os¬ 
borne et al. 1982,1984; Alexander et al. 1984; Hilgemeier et al. 1988) and three 
that are based on the late decay of ^Be (Osborne et al. 1982,1984; Robertson et 
al. 1983; Volk et al. 1983). There is a systematic discrepancy of more than 3a 
between these two data sets whose origin is not clear. When theoretical mod¬ 
els are used to extrapolate the direct measurement to low energies they yield 
a weighted mean 834 ( 0 ) = 0.507 ± 0.016 MeV • b, while the activation measure¬ 
ments yield 834 ( 0 ) = 0.579 ± 0.024 MeV • b. The weighted mean of all experi¬ 
ments, 834 ( 0 ) = 0.53 ± 0.05 (Adelberger et al. 1998) was used in BSP98. How¬ 
ever, the radiative capture ^He He Be -b 7 does not take place well out¬ 
side the range of the nuclear forces jE b ~ 3.6 /m < R(^He) -b R(^He)). 

Thus, the ability of the theoretical models to predict correctly the low energy 
dependence of 834 (E) is questionable. Moreover, the energy dependence of the 
cross section which they predict disagrees with that observed in the measure¬ 
ments of Krawinkel et al. 1982 who used gas targets, have low statistical errors 
and extend to the lowest energies. If one uses a polynomial fit of S 34 (E) to these 
data, and than uses it to extrapolate the other measurements to low energies, one 
obtains from the direct measurements a weighted mean 834 ( 0 ) = 0.45 MeV • b 
which was used in D896. 

6.3 Diffusion 

Diffusion, caused by density, temperature, pressure, chemical composition and 
gravitational potential gradients plays an important role in the sun since it 
modifies the local chemical composition in the sun. The agreement between 
helioseismology and the 88 M is improved when diffusion is included (see e.g.. 
Fig. 10), but inclusion increases the discrepancy between the 88 M and the solar 
neutrino observations. The relative changes in 88 M predictions due to diffusion 
of all elements are summarized in Table IV. While BP95 and BSP98 found a 
rather large increase in the predicted ^Be, ®B, and solar neutrino 

fluxes; 14%, 36%, 52%, 58%, and 61% which result in 36%, 33%, 9% increases 
in their predicted rates in 8 uper-Kamiokande, Homestake, and in GALLEX and 
SAGE, respectively, DS96 found only a moderate increase due to diffusion, 4%, 
10%, 23%, 24% and 25%, respectively, in the above fluxes, which result in 10%, 
10% and 2% increase in the predicted rates in Kamiokande, Homestake, and 
in GALLEX and SAGE, respectively. Although the two diffusion calculations 
assumed a different initial solar chemical composition (see below) and BP95 
approximated the diffusion of all elements heavier than ^He by that of fully ion- 
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ized iron (the DS calculations followed the diffusion of each element separately 
and uses diffusion coefficients calculated for the local ionization state of each 
element in the sun as obtained from solving the local Saha equations), these 
cannot fully explain the above large differences. Independent diffusion calcu¬ 
lations by Richard et al. (1996) obtained similar results to those obtained in 
DS96 as can be seen from Table IV (we interpolated the results from the two 
models of Richard et al. (1996) to the initial chemical composition assumed in 
DS96). Note that internal magnetic fields can suppress diffusion significantly. 

7 New Neutrino Physics? 

Standard solar models, like DS96, perhaps can explain the results reported by 
Kamiokande and Super-Kamiokande. However, if the neutrino absorption cross 
sections assumed by the radiochemical experiments are correct, then present 
standard solar models cannot explain the absence of the expected contributions 
of the ^Be and CNO solar neutrinos to the ^^Ar production rate in ^^Cl and 
to the ^^Ge production in ^^Ga. Gonsequently, many authors have claimed 
that the the solar neutrino observations imply neutrino properties beyond the 
minimal standard electroweak model (e.g., Bahcall and Bethe 1991). 

Neutrino Magnetic Moments. Some authors have suggested that neutrinos 
may have anomalous magnetic moments large enough so that the solar internal 
magnetic field in the sun can flip the neutrino helicity and convert part of the 
left handed weakly interacting solar neutrinos into noninteracting right handed 
neutrinos (e.g., Okun et al. 1986, Lim and Marciano 1988, Akhmedov 1988). 
However, the high statistics solar neutrino measurements of SK do not show 
the time variation of the solar neutrino flux which is predicted by a magnetic 
helicity flip interpretation of the solar neutrino anomaly. 

Neutrino Oscillations: Mikheyev and Smirnov (1985) have discovered that 
neutrino oscillations in matter (Wolfenstein 1978, 1979) can lead to resonant 
conversion of neutrino flavor (the MSW effect) in the sun and explain the so¬ 
lar neutrino observations quite neatly. It requires only a natural extension of 
the minimal standard electroweak theory and it is based on a simple quan¬ 
tum mechanical effect. Many authors have carried out extensive calculations to 
determine the neutrino mixing parameters which can bridge between the pre¬ 
dictions of the standard solar models and the solar neutrino observations. The 
neutrino mixing parameters can also be deduced analytically directly from the 
solar neutrino observations (e.g., Dar and Nussinov 1991; Dar 1993): If the v,. 
is mixed with the (or the with a vacuum mixing angle 0 <C I and a mass 
difference Am^ = then solar i/g’s which are produced in the sun can 

flip their flavor on their way out of the sun if they encounter electron density 

Am‘^c^cos20 Aw?c^ 

ne =- 1= -«—^-. (17) 

2V2GfE,^ 2V2GfE,^ 
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The probability for the resonant flavor flip (the MSW effect) is given approx¬ 
imately by (e.g., Haxton 1986; Parke 1986; Dar et al. 1987) P{ve ^ ~ 

1 — where 


t:H Am‘^c^sin^29 -nHAm^c^O^ 

Ahc ■ cos29 Tic ’ 

with H = —ne/{dne/dr) being the density scale-height at the resonance. 

Strong suppression of the contribution from the pep, ^Be and CNO solar 
neutrinos to the Cl experiment requires a complete flavor change of these neu¬ 
trinos in the sun, i.e., that they encounter a resonant density with e ;g> Ei,^ ~ 1 
MeV. But the results from GALLEX and SAGE suggest that the lower energy 
pp solar neutrinos evade such a flavor flip. This is possible if the resonance con¬ 
dition (17) is satisfied for the pep, ^Be and CNO solar neutrinos (Aj^^ > 0.862 
MeV) but not for the pp neutrinos (Aj^^ < 0.420 MeV). Since the maximal (cen¬ 
tral) electron density in the sun is Ue « IO^Na the last condition reads (e.g., 
Dar 1993) 

0.5 X 10“®eV^ < Arn^c'^ < 1 x lO'^eV^ (19) 

Eor a given Arri^, eq. (18) can be used to adjust the mixing angle to produce 
the required suppression of the SSM ®B solar t'e flux: 

Am^sin^29 « 3 x 10“®eV^c“"^, Am^sm^20 « 1 x 10“®eV^c“^, (20) 

for respectively, the BSP98 and the DS96 standard solar models. 

8 Are The ^Be Solar Neutrinos Missing ? 

Electron capture by ^Be into the ground state of ^Li produces 862 keV neutri¬ 
nos. The threshold energy for neutrino absorption by ^"^Cl is 814 keV. Thus, 
absorption of ^Be neutrinos by ^^Cl produces 48 keV electrons. The maximum 
energy of the pp solar neutrinos is 420 keV. The threshold energy for neutrino 
absorption in ^^Ga (3/2 “) is 233 keV into the ground state (l/2“) and 408 into 
its first excited state (5/2“). The produced electrons have therefore energies 
below 187 and 12 keV, respectively. If the theoretical cross sections for neu¬ 
trino absorption near threshold overestimate their true values significantly then 
the predicted rates will significantly overestimate the expected signals in the 
Chlorine and Gallium experiments. An indication that final state interactions 
effects are not completely understood is provided by Tritium /3-decay. Although 
final state interactions in Tritium /3-decay have been studied extensively, they 
do not explain well the end-point /3-decay spectrum {Ee ~ 18.6 keV). In all 
recent measurements, the measured spectrum yields a negative value for the 
fitted squared mass of the electron neutrino. Final state interactions effects 
(screening of the nuclear charge by atomic electrons, exchange effects, radiative 
corrections, nuclear recoil against the electronic cloud, etc) in neutrino captures 
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near threshold in and Ga may be much larger because their Z values are 
much larger and because the de Broglie wave lengths of the produced electrons 
are comparable to the Bohr radii of the atomic K shells in Cl and Ga. If final 
state interactions reduce considerably the near threshold absorption cross sec¬ 
tions of pp neutrinos in ^^Ga (making room for the expected contribution of 
^Be solar neutrinos in Gallium) and of ^Be neutrinos in ^^Cl, perhaps they can 
make the solar neutrino observations of Kamiokande and the Homestake exper¬ 
iment compatible. Such an explanation of the solar neutrino problem implies 
that experiments such as BOREXINO and HELLAZ will observe the full ^Be 
solar neutrino flux. 


9 Astrophysical Solutions To The SNP 


Even if the ^Be solar neutrino flux is strongly suppressed, it does not eliminate 
standard physics solutions to the solar neutrino problem: The ratio between the 
fluxes of ^Be and ®B solar neutrinos is given by 

^ _ <(>t/ 0 (^Be) _ f UeUT < av >e 7 d^r 
f UpUr < av >p 7 (Pr 


Because of the decreasing temperature and Be7 abundance as function of dis¬ 
tance from the center of the sun on the one hand, and the ~ increase in 
radial mass on the other, the production of ^Be and ®B solar neutrinos in the 
SSM peaks around an effective radius, Veff ~ O.O64i?0 (re// is approximately 
the radius within which 50% of the flux is produced) . The SSM also predicts 
a ratio of electron to proton densities near the center of the sun, rie/rip ~ 2 , 
consistent with helioseismology observations. Gonsequently, the SSMs predict 


2 < av >e 7 
< av >p 7 


1.27 X 




( 22 ) 


where F 17 is the screening correction to the p-capture rate by ^Be, T 7 is the 
temperature in lO^RT at the effective radius and S 17 is in eV barn units. The 
SSMs yield T 7 (re//) « 1.45. Using S'i 7 ( 0 ) = 17 eU b and (() 0 (®B) = 2.44 x 
10® cm~^ s~^ as observed by Super-Kamiokande, one can reproduce the SSM 
prediction (e.g., Dar and Shaviv 1996) 

<^^g(7Be) = R(|)^,QCB) « 3.7 x 10® cm"® (23) 


Astrophysical solutions of the solar neutrino problem aim towards suppressing 
the value of R. Three alternatives are currently investigated: 

Plasma Physics Effects: The effects of the surrounding plasma on nuclear 
reaction rates in dense stellar plasmas, and in particular on proton and electron 
capture by ^Be in the sun are known only approximately. In order to explain 
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the deficit of "^Be solar neutrinos, without much affecting the SSM, plasma 
screening effects must reduce considerably the electron/proton capture ratio by 
^Be, relative to the predictions of the weak screening theory (Salpeter and Van 
Horne 1969). The screening enhancement of bare nuclear cross sections is not 
well understood even in laboratory measurements with gas targets. Also the 
applicability of the weak screening theory to the dense plasma in the solar core 
is questionable. Moreover, correlations and fluctuations, which are neglected 
in the weak screening theory can affect strongly the screening enhancement of 
nuclear reaction rates in the solar core. This possibility is currently studied, e.g., 
by Shaviv and Shaviv (1998) using numerical methods. Because of accidental 
cancellations the weak screening corrections to the rates of all nuclear reactions 
do not change the predicted ®B solar neutrino flux, but perhaps a more exact 
treatment of screening may change R considerably. 

Because the sub-Coulomb nuclear reactions in the core of the sun take place 
mainly between nuclei with kinetic energies much larger than their mean ki¬ 
netic energies, their rates are very sensitive to the high energy tail of their 
velocity distribution in the sun. Diffusion, radiative flows, energetic nuclear 
products, internal fluctating (equipartition ?) electric and magnetic fields and 
other collective effects may change the assumed Maxwell-Boltzmann tails of the 
energy distribution of the energetic particles in the core of the sun. This may 
shift the position of the Gamow peaks for the nuclear reaction rates and change 
considerably the ratios between nuclear reaction rates in the sun which have 
very different temperature dependence (see, e.g., Kaniadakis et al. 1998). 

In principle, collective plasma physics effects, such as very strong magnetic or 
electric fields near the center of the sun, may also polarize the plasma electrons, 
and affect the branching ratios of electron capture by ^Be (spin 3/2“) into the 
ground state (spin 3/2“, = 0.863 MeV, BR=90%) and the excited state 

(spin l/2“, = 0.381 MeV, BR=10%) of ^Li. Since solar neutrinos with 

Eue = 0.381 MeV are below the threshold (0.81 MeV) for capture in ^^Cl and 
have a capture cross section in ^^Ga that is smaller by about a factor of 6 
relative to solar neutrinos with E,^^ = 0.863 MeV, therefore a large suppression 
in the branching ratio to the ground state can produce large suppressions of the 
^Be solar neutrino signals in ^^Cl and in ^^Ga. However, such an explanation 
requires anomalously large fields near the center of the sun. 

Temporal and Spatial Variations in T: Davis (1996) has been claiming 
persistently that the solar neutrino flux measured by him and his collabora¬ 
tors in the ^^Gl radiochemical experiment is varying with time. Because of the 
possibility that neutrinos may have anomalous magnetic moments, much larger 
than those predicted by minimal extensions of the standard electroweak model, 
which can solve the solar neutrino problem, attention has been focused on an¬ 
ticorrelation between the solar magnetic activity (the 11 year cycle) and the 
vq flux (see, e.g., Davis 1996). Also a day-night effect (e.g., Cribier et al 1986; 
Dar and Mann 1987) due to resonant conversion of the lepton flavor of solar 
neutrinos which cross Earth at night before reaching the solar neutrino detector 
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was not found by Kamiokande. However, the basic general question whether 
the solar neutrino flux varies on a short time scale, has not been fully answered, 
mainly because of the limited statistics of the first generation of solar neutrino 
experiments. The SSM predict no significant variation of the solar neutrino flux 
on time scales shorter than millions of years. However, the sun has a differential 
rotation. It rotates once in ~ 25 days near the equator, and in ~ 33 days near 
the poles. Moreover, the observed surface rotation rates of young solar-type 
stars are up to 50 times that of the sun. It suggest that the sun has been losing 
angular momentum over its lifetime. The overall spin-down of a sun-like star 
by mass loss and electromagnetic radiation is difficult to estimate from stellar 
evolution theory, because it depends on delicate balance between circulations 
and instabilities that tend to mix the interior and magnetic fields that retard or 
modify such processes. It is quite possible that the differential rotation extends 
deep into the core of the sun and causes there spatial and temporal variations in 
the solar properties due to circulation, turbulences and mixing. Since R is very 
sensitive to the temperature, even small variations in temperature can affect R 
significantly without affecting significantly the pp solar neutrino flux (the ^Be 
and ®B solar neutrinos will come mainly from temperature peaks, while the pp 
neutrinos will reflect more the average temperature). If the solar neutrino flux is 
time dependent, then cross correlation analysis of the various data sets from the 
Homestake, Kamiokande, GALLEX, SAGE and Superkamiokande may reveal 
such unexpected correlations: If arbitrary time lags are added to the different 
solar neutrino experiments, the cross correlation is maximal when these time 
lags vanish. Moreover, a power spectrum analysis of the signals may show peaks, 
if the time variation is periodic. In particular, Super-Kamiokande with its high 
statistics should examine whether data from different fiducial volumes are cross 
correlated in time. Relevant information on time variability in the solar core 
may come soon also from SOHO and GONG. 

Mixing of ^He: The SSM ^He equilibrium abundance increases sharply with 
radius. Cummings and Haxton (1996) have recently suggested that the ^Be so¬ 
lar neutrino problem could be circumvented in models where ^He is transported 
into the core in a mixing pattern involving rapid filamental flow downward. We 
note that if this mixing produces hot spots (due to enhanced energy release) 
they can increase the effective temperature for p capture by ^Be in a cool en¬ 
vironment, reducing R while keeping the ®B solar neutrino flux at the observed 
level. Perhaps, helioseismology will be able to test that. 

Cummings and Haxton (1996) also noted that such mixing will have other as- 
trophysical consequences. Eor example, galactic evolution models predict ^He 
abundances in the presolar nebula and in the present interstellar medium (ISM) 
that are substantially (i.e., a factor of five or more) in excess of the observation- 
ally inferred values. This enrichment of the ISM is driven by low-mass stars in 
the red giant phase, when the convective envelope reaches a sufficient depth to 
mix the ^He peak, established during the main sequence, over the outer portions 
of the star. The ^He is then carried into the ISM by the red giant wind. The 
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core mixing lowers the main sequence ^He abundance at large r. 


10 Conclusions 

The solar neutrino problem does not provide conclusive evidence for neutrino 
properties beyond the standard electroweak model. The solar neutrino problem 
may be an astrophysical problem. The deviations of the experimental results 
from those predicted by the standard solar models may reflect the approximate 
nature of our knowledge of nuclear reaction rates and radiation transport in 
dense stellar plasmas and the approximate nature of the solar models which 
neglect angular momentum effects, differential rotation, magnetic field, angular 
momentum loss and mass loss during evolution and do not explain yet, e.g., so¬ 
lar activity and the surface depletion of lithium, relative to its meteoritic value 
(which may or may not be relevant to the solar neutrino problem). Improve¬ 
ments of the standard solar model should continue. In particular, dense plasma 
effects (screening, correlations, fluctuations and deviations from Maxwell- Boltz¬ 
mann distributions) on nuclear reaction rates and radiative opacities, which 
are not well understood, may affect the SSM predictions and should be fur¬ 
ther studied, both theoretically and experimentally. Relevant information may 
be obtained from studies of thermonuclear plasmas in inertial confinement ex¬ 
periments. Useful information may also be obtained from improved data on 
screening effects in low energy nuclear cross sections of ions, atomic beams and 
molecular beams incident on a variety of gas, solid and plasma targets. 

Better knowledge of low energy nuclear cross sections is still needed. Im¬ 
proved measurement of the low energy nuclear cross sections for the radiative 
captures p Be B -|- 7 and ^He +‘^ He Be -I- 7 by photodissociation of 
®B and ^Be in the Coulomb field of heavy nuclei can help determine whether 
there is a ®B solar neutrino problem. 

Neutrino oscillations, and in particular the MSW effect, may be the correct 
solution to the solar neutrino problem. But, only future experiments, such as 
SNO, BOREXINO and HELLAZ, will be able to supply a definite proof that 
Nature has made use of this beautiful effect. 
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Table la: Comparison between the solar neutrino fluxes predicted by the SSM 
of BP95, BSP98, BTM98 and DS96, and measured by the solar neutrino exper¬ 
iments. 


V Flux 

BP95 

BSP98 

BSP98 

DS96 

Experiment 

(j)v{pp) 

5.91 

5.94 


6.10 


4)y{pep) [10® 

1.39 

1.39 


1.43 


(j)SCBe) [10®CTO“^s“^] 

5.18 

4.80 


3.71 


4>v^B) [10®CTO“^S“^] 

6.48 

5.15 

4.82 

2.49 

2.44 ±0.11 

[10® 

6.4 

6.05 


3.82 


[10® 

5.15 

5.32 


3.74 


[10®CTO“^s“^] 

6.48 

6.33 


4.53 


Rates 






E(0a)ci [SNU] 

9.3 ± 1.4 

7.7 ± 1.2 

7.18 

4.1 ± 1.2 

2.56 ±0.25 

E(0a)Ga [SNU] 

137 ±8 

129 ±8 

127 ±8 

115 ±6 

76.4 ± 8 

n^a)Ga [^iVt/] 

137 ±8 

129 ±8 

127 ±8 

115 ±6 

70 ±8 


Table Ib Characteristics of the BP95, BTM98, and DS96 Solar Models in Table 
la (c=center; s=surface; bc=base of convective zone; N = log([N]/[H]) -|- 12). 


Parameter 

BP95 

BTM98 

DS96 

T, [lO^iF] 

1.584 

1.567 

1.561 

Pc [g cm”®] 

156.2 

151.9 

155.4 


0.3333 

0.3442 

0.3424 

Me 

0.6456 

0.635 

0.6380 

Zc 

0.0211 

0.02084 

0.01940 

Rconv [R/Rq] 

0.712 

0.715 

0.7130 

The [10®K] 

2.20 

2.172 

2.105 

X, 

0.7351 

0.739 

0.7512 


0.2470 

0.243 

0.2308 


0.01798 

0.0181 

0.0170 

NsC^C) 

8.55 

8.55 

8.55 

Nsi^*N) 

7.97 

7.97 

7.97 

NsC^O) 

8.87 

8.87 

8.87 

Ns^Ne) 

8.08 

8.08 

8.08 

Teff [K] 

5800 

5800 

5803 
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Table II: Comparison between the SSM of Bahcall and Pinsonneult (1995) and 
of Dar and Shaviv (1996). 


Parameter 

BSP98 

DS96 

Mq 

1.9899 X 10^3 g 

1.9899 X 10®3 

Lq 

3.844 X 10^^ erg 

3.844 X 10^3 e 

Rq 

6.9599 X 10^° cm 

6.9599 X lO^o 

Iq 

4.566 X 10® y 

4.57 X 10® y 

Rotation 

Not Included 

Not Included 

Magnetic Field 

Not Included 

Not Included 

Mass Loss 

Not Included 

Not Included 

Angular Momentum Loss 

Not Included 

Not Included 

Premain Sequence Evolution 

Initial Abundances : 

Not Included 

Included 

^He 

Adjusted 

Adjusted 

C,N,0,Ne 

Adjusted 

Adjusted 

All Other Elements 

Photospheric Abundances : 

Adjusted 

Meteoritic 

^He 

Predicted 

Predicted 

C,N,0,Ne 

Photospheric 

Photospheric 

All Other Elements 

Meteoritic 

Predicted 

Radiative Opacities 

OPAL 1996 

OPAL 1996 

Equation of State 

Straniero 1996? 

DS 1996 

Partial Ionization Effects 
Diffusion of Elements : 

Not Included 

Included 

H, ^He 

Included 

Included 

Heavier Elements 

Approximated by Fe 

All Included 

Partial Ionization Effects 
Nuclear Reaction Rates : 

Not Included 

Included 

Sii(O) eV • b 

4.00 X 10"^® 

4.07 X 10-39 

S33(0) MeV-b 

5.3x 

5.6 

S34(0) keV • b 

0.53 

0.45 

Si7(0) eV • b 

19 

17 

Screening Effects 

Included 

Included 

Nuclear Equilibrium 

Imposed 

Not Assumed 
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Table IV: Fractional change in the predicted vq fluxes and counting rates in the 
uq experiments due to the inclusion of element diffusion in the SSM calculations 
of Bahcall and Pinsonneault (BP95), Brun, Turck-Chieze and Morel (BTM98), 
Dar and Shaviv (DS96) and Richard, Vauclair, Charbonnel and Dziembowski 
(RVCD96). The results of models 1 and 2 of RVCD96 were extrapolated to the 
initial solar composition which was used in DS96. 



BP95 

BTM98 

DS96 

RVCD96 

PP 

- 1.7% 



- 0.3% 

- 0.8% 

pep 

- 2.8% 



- 0.3% 

- 0.4% 

^Be 

+13.7% 



+4.2% 

+ 6.5% 

®B 

+36.5% 

+32 

% 

+11.2% 

+ 10.7% 

13n 

+51.8% 



+22.7% 

+ 19.8% 

150 

+58.0% 



+24.0% 

+20.8% 

ITp 

+61.2% 



+24.9% 

+21.8% 

Rates 






H20 

+36.5% 

+32 

% 

+11.2% 

+ 13.3% 

Cl 

+32.9% 

+27 

% 

+ 9.5% 

+ 12.3% 

Ga 

+ 8.7% 

+ 7 

% 

+ 2.6% 

+ 3.7% 


26 



Figure Captions 

Fig. 1. The principal branches of the pp cycle and the CNO bi-cycle. 

Fig. 2. The solar neutrino capture rate in ^^Cl as measured in the Homestake 
experiment runs nos. 18 tO 133. The dashed line shows the average value. 

Fig. 3. The ®B solar neutrino flux as function of time as measured by 
Kamiokande between 1986-1996. 

Fig. 4. The cosine of the angle between the electron direction and the radius 
vector from the sun in Super-Kamiokande. The solid line shows the best fit for 
a ®B solar neutrino flux. 

Fig. 5. The time variation of the ®B solar neutrino flux as measured by 
Super-Kamiokande from June 96 to June 97. The solid line shows the expected 
variation of the flux due to the eccentricity of the Earth’s orbit around the sun. 

Fig. 6. The ratio between the ®B solar neutrino flux observed by Super- 
Kamiokande and the flux predicted by the SSM of BP95 as function of zenith 
angle of the sun. 

Fig. 7. The ratio between the observed number of electrons scattered by solar 
neutrinos in Super-Kamiokande and their expected number in the SSM of BP95 
as function of electron recoil energy. 

Fig. 10 The ^^Ge production rate in ^^Ga by solar neutrinos as measured by 
GALLEX and SAGE between 1990 and 1997. 

Fig. 9. The astrophysical S(E) factor for the reaction ^He -|-^ He He -I- 2p 
as measured by various low energy experiments. The dotted line is a best 
polynomial fit (solid line) to the data with maximal screening enhancement 
(C/e = 240 eV). 

Fig. 10. The relative difference beteen the speed of sound squared as in¬ 
ferred from the helioseismological measurements of GOLF and LOWL on board 
SOHO and that calculated in BTM98 from the SSM with diffusion (solid line) 
and without diffusion (dotted line). 
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